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We study J/tp production at RHIC and LHC energies with both initial production and regener- 
ation. We solve the coupled set of transport equation for the J/tp distribution in phase space and 
the hydrodynamic equation for evolution of quark-gluon plasma. At RHIC, continuous regeneration 
is crucial for the J/tp momentum distribution while the elliptic flow is still dominated by initial 
production. At LHC energy, almost all the initially created J/tps are dissociated in the medium and 
regeneration dominates the J/lj) properties. 
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PACS numbers: 25.75.-q, 12.38.Mh, 24.85.+p 



The goal of high-energy nuclear collisions is to iden- 
tify and study the equation of state of the Quark-Gluon 
Plasma (QGP) which is believed to exist at the early 
stage of our universe In ^/snn— 200 GeV Au + 
Au collisions, the observations of jet-quenching Q and 
collective flow have demonstrated the formation of hot 
and dense matter with partonic collectivity [1,01 ■ Local 
thcrmalization of the system created in heavy ion colli- 
sions is yet to be tested. Heavy flavors including charm 
and bottom quarks are powerful tools [f| for studying 
the early collision dynamics because their masses are 
much larger than possible excitation, temperature for 
example, of the system created in collisions at RHIC. 
Recent results have shown that heavy flavors are all pro- 
duced in the initial collisions Therefore, study- 
ing open charm and charmonium production will yield 
important information on the properties of QGP. The 
J/tp is a particularly sensitive probe of the early stages 
because its survial probability depends on the environ- 
ment. Lattice gauge theory calculations have indicated 
that J / tps do exist above the critical temperature [1, Q . 

The J/tp suppression was first proposed as a direct 
signature to identify the QGP formation twenty years 
ago • Besides the normal suppression induced by nu- 
clear absorption, the J ftps initially produced by hard 
processes are anomalously suppressed [III EH EH EH 
EH EH E3 by interactions in the hot medium. While 
charm quark production at the SPS is expected to be 
small, there are more than 10 cc pairs produced in a 
central Au + Au collision at RHIC, and is probably 
more than 200 pairs at the LHC [Tg] - These uncorre- 
cted charm quark pairs in the QGP can be recombined 
to form J /tps. Obviously, regeneration will enhance the 
J/tp yield and alter its momentum spectra. Recently, 
the regeneration approach for J/tp production at RHIC 
has been widely discussed with different models, such as 
thermal production on the hadronizaton hypersurface 
according to statistic law EH HH HH , the coalescence 
mechanism [22j and the kinetic model [H, [24[ which 
considers continuous J/tp regeneration in a QGP. 

The medium created in high-energy nuclear collisions 
evolves dynamically. In order to extract information 
about the medium by analyzing the J/tp distributions, 
both the hot and dense medium and the J/tp produc- 
tion processes must be treated dynamically. In this 
Letter, we treat continuous regeneration of the J/tp in 
a QGP self-consistently, including hydrodynamic evo- 
lution of the QGP itself. Both regenerated and initially 



produced J /tps are simultaneously taken into account. 
Including both processes is important for RHIC since, 
unlike the collisions at the LHC, not all the initially 
produced J /tps are destroyed. 

The massive J/tps are unlikely fully thermalized 
with the medium. Thus their phase space distribution 
should be governed by transport rather than the kinetic 
equation of Refs. [H, [24[. To comprehensively treat 
the J/tp distribution, we determine the J/tp transport 
equation including both initial production and anoma- 
lous suppression as well as regeneration. The transport 
equation is then solved together with the hydrodynamic 
equation which characterizes the space-time evolution 
of the QGP. In order to compare with the experimen- 
tal measurements, we calculate the average transverse 
momentum squared, (p 2 ), elliptic flow, t>2, and nuclear 
modification factor, Raa, for J/tp at both RHIC and 
LHC. 

In p + p collisions, about 30-40% of final state J/tps 
are from the feed-down [25| of tp' and \c ■ In order 
to simplify the numerical calculation, we neglect the tp' 
contribution and assume 40% of the final state J ftp's 
come from \c ">'• Since the $ (=J/tp,Xc) is heavy, 
we use a classical Boltzmann-type transport equation 
to describe its evolution. The distribution function 
/*(Pti x t? T |b) m t ne central rapidity region and in the 
transverse phase space (pt,x t ) at fixed impact param- 
eter b is controlled by the equation 



(1) 



The second term on the 
free-streaming of \& with 



left-hand side arises from 
transverse velocity = 



Pt/vP* + TO *i which leads to the "leakage" effect and 
is needed to explain the averaged transverse momen- 
tum squared at the SPS [15j . The anomalous suppres- 
sion and regeneration mechanisms are reflected in the 
loss term and gain term ply, respectively. 

Suppose the medium locally equilibrates at time To, 
when nuclear absorption of the initially produced J/tps 
has ceased. Absorption effects can be included in 
the initial distribution [IHEH], /*(pt, x*, ro|b), of the 
transport equation. To determine the initial condition 
at RHIC, we take r = 0.6 fm, the initial J/tp produc- 
tion cross section a pp = 2.61 ± 0.20 ub 12611 and nuclear 
absorption cross section a a bs = 3 mb [27] . We also as- 
sume that a gN = 0.076 GeV 2 /fm and (pt) pp = 4.31 

GeV 2 (2(| to describe the initial p t broadening due to 
multiple gluon scattering. 
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When the loss and gain terms a* and are known, 
the transport equation ([T]) can be solved analytically 
with the result 



/*(pt,x t ,r|b) = /*(p t ,x t -v*(r-To),r |b) 



(2) 



x e ~ /to dT'a*(pt,xt-v*(T-T'),r'|b) 



/ rfr'/3*(pi,x t - v^(r - r'),T'|b) 

•/ Tn 



x e 



- dT"a*(p t ,x t -v*(T-T"),r"|b) 



The first and second terms on the right-hand side in- 
dicate the contributions from the initial production 
and continuous regeneration, respectively. Both suf- 
fer anomalous suppression. The coordinate shift x f — > 
Xi — v$At reflects the leakage effect during the time 
period At. 

We determine now the J/tjj suppression and regener- 
ation in a QGP. We consider only the gluon dissociation 
process, g + ty — > c + c, for the loss term [lfij]. Its inverse 
process is the gain term. Then a$ and are 



a*(p t ,xj,r|b) 



/3*(p t ,xt,r|b) 



1 



d 3 Pg 



2E^ J (2n) 3 2E g 
x e(T(x t ,r|b)-T c ), 
1 /' ri 3 p ff d 3 p c 



Wgl(a)/g(Ps,x t , 



d 3 t 



7 (27r) 3 2£ g (2 7 r) 3 2 J B c (27r) 3 2 J B e 
x W£*(«)/ C (p c , x t) r|b)/ g (p gj x t , r|b) 



x(27r) 4 S^(p + p g - Pc 
x6(T (x t ,r|b)-T c ), 



(3) 



where E g , Ey, E c and -Eg are the gluon, charmonium, 
c and c energies. The step function Q ensures that 
anomalous suppression and regeneration occur only in 
the QGP phase. The gluon thermal distribution is f g . 
W^p( s ) is transition probability of the gluon dis- 
sociation as a function of s = (p + p g ) 2 calculated in 
pQCD [28|, and Wf& is c and c recombination tran- 
sition probability. Note that Wff can be obtained 
from using detailed balance. In numerical calcu- 
lations, we take m c = 1.87 GeV, including in-medium 
effects HE H|, m J/t = 3.1 GeV, m Xc = 3.51 GeV and 
critical temperature T c = 165 MeV. 

Unlike gluons that are constituents of a QGP, charm 
quarks are heavy and may not be thermalized in the 
QGP. In principle, a self-consistent treatment of charm 
quark motion in the QGP should be described by a 
transport equation similar to Eq. |T]). For simplicity, 
we consider the c and c distribution functions, f c and 
fc, in two extreme scenarios. In the weak interaction 
limit, the charm quarks in QGP are assumed to keep 
their original momentum distribution, g c (q), calculated 
in pQCD [24[ and their initial space distribution deter- 
mined by nuclear geometry, 

/ C)g (q,x t |b) = af p T A (x t )T B (x t - b) 5c (q), (4) 

where a pp = 622 ± 57 /xb [1| is the charm production 
cross section at RHIC. In the opposite limit, when the 
c and c are strongly correlated to the medium, they 
are thermalized and distributed statistically. Both the 
pQCD and thermal distributions are normalized to the 
initial charm quark number. 



Since the hadronic phase occurs later in the evolution 
of heavy ion collisions when the density of the system 
is lower compared to the early hot and dense period, 
we have neglected the hadronic dissociation. The sup- 
pression and regeneration region controlled by the step 
function in the loss and gain terms and the gluon and 
thermal charm quark distributions are determined by 
the QGP evolution. With the Hubble-like longitudi- 
nal expansion and boost invariant initial condition, the 
transverse hydrodynamical equations [16] can be solved 
numerically to determine the evolution of temperature 
and transverse fluid velocity. 

The final state J/ip distributions at fixed central- 
ity can be found by integrating the distribution func- 
tion, Eq. ([5]), over p t and x< up to r — > oo. We first 
examine the mid-rapidity nuclear modification factor 
Raa- The numerical result, as a function of partici- 
pant number, N part , is compared to the mid-rapidity 
RHIC data in Fig. [TJ Plots (a) and (b) correspond 
to the pQCD and thermal charm quark distributions, 
respectively. The bands are theoretical results includ- 
ing the uncertainty in the initial charm quark and VE' 
production cross sections and a^ p . The solid and 
dot-dashed curves indicate calculations with only initial 
production (without regeneration, = 0) and only 
regeneration (without initial production, a a bs = oo). 
Almost all the initial 40% of J ftps from \c decay is 
lost in semi-central collisions [3l|. In central collisions, 
only the directly produced J/ips suffer anomalous sup- 
pression. When nuclear absorption effect is reduced by 
decreasing the initial time To or the absorption cross 
section a a bs, the initial J/tp yield is less suppressed. 
Due to the strong regeneration in central collisions, the 
full J ftp yield is no longer a monotonically decreasing 
function of centrality. The data seem to show a flat 
region at 50 < N part < 150. This feature can not be 
reproduced in the present model calculations. 
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FIG. 1: The nuclear modification factor Raa as a function 
of participant number N par t in pQCD (a) and thermal (b) 
charm quark distributions. The solid and dot-dashed curves 
are the calculations with only initial production and only 
regeneration. The bands are the full result, including the 
uncertainty in initial charm quark and ^ production. The 
data are from PHENIX 

The momentum spectra are expected to be more sen- 
sitive to the production mechanism than the integrated 
yields. We show the J/ip average squared transverse 
momentum, (p 2 ), as a function of centrality in Fig. [2] 
It is well known that the multiple gluon scattering in 
the initial state leads to transverse momentum broad- 
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ening. The leakage effect due to the anomalous sup- 
pression also results in p t broadening in central colli- 
sions 15]. These two broadening effects on the initial 
J ftps are reflected in the solid curves in Fig. [2J For 
the regenerated J ftps, the charm quarks in the pQCD 
scenario have no scattering in the initial state and in 
the QGP, and they satisfy the statistic distribution in 
the thermal scenario. The Opf) of the regenerated J/ips 
in these two limits looks independent of centrality, see 
dot-dashed curves in Fig. [5J Since the high momentum 
charm quarks lose energy during thermalization, the 
average squared momentum (~ 2 GcV 2 ) in the ther- 
mal scenario is smaller than the one (~ 2.5 GeV 2 ) 
in the pQCD scenario. The collective flow develops 
with time, the (p 2 ) for the sudden-produced J/ips on 
the hadronization hypersurface [HI, [2(1 H2, [H| of QGP 
is larger than our result for continuously regenerated 
J ftps in the whole volume of QGP. Because the initial 
production decreases with centrality, due to absorption 
and anomalous suppression, and regeneration increases 
with centrality, shown in Fig. [TJ the (p 2 ) starts at the 
initial result and approaches the regeneration result as 
centrality increases. 
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FIG. 2: The J/ip averaged transverse momentum square 
as a function of the number of binary collision, N C oll- The 
calculations are the same as in Fig. [TJ The PHENIX [3^] 
data are shown. 

The J/ip elliptic flow, is a useful tool for study- 
ing charm production mechanism in nuclear collisions. 
When no regeneration is considered, the leakage effect 
gives the lower limit of J/ip v 2 [la ], the solid curve in 
Fig. [3J In the pQCD scenario, charm quarks do not 
interact with the medium and the regenerated J/ips 
do not carry any collective property of the system, so 
we only discuss v 2 in the thermal scenario. If the re- 
generation occurs only at hadronization of the QGP, 
the well-developed azimuthal anisotropic flow leads to 
a large J/ip v 2 , see discussions in Ref. [22j |. However, 
for the continuous regeneration in the QGP volume the 
average elliptic flow becomes much smaller since the 
asymmetric flow is very small in the early stage. At 
impact parameter b = 7.8 fm, the continuously regen- 
erated J/ip v 2 , the dot-dashed line in Fig. [3J is only 
one half of the suddenly regenerated J/ip v 2 in the co- 
alescence model [22j]. To check the difference between 
the continuous and sudden regeneration mechanisms, 
we force the regeneration to occur in a thin spherical 
shell by adjusting the regeneration temperature region 
to T c < T < T c + ST. Indeed, the regenerated v 2 ap- 
proaches the value of Ref. [12] for 5T — > 0. Since the 



J/ip yield in semi-central collisions is still dominated 
by the initial production, see Fig. [TJ there is no siz- 
able difference in v 2 between the total (dashed line) 
and initial-produced J/ips (solid line), see Fig.[3J 

What is the case for nuclear collisions at the LHC? 
Due to the extremely high center of mass energy, the 
QGP formed at LHC energy will have much higher tem- 
peratures, longer lifetimes and larger sizes. Most of 
the initial J/ips will be suppressed by stronger gluon 
dissociation in QGP. On the other hand, more char- 
monia will be created by the regeneration. The ratio 
N reg /N ml of regenerated to initial-produced J/ips is 
shown as a function of centrality in Fig. [4] While at 
RHIC, the ratio can reach unity only in the most cen- 
tral collisions, it is much larger than one at almost any 
centrality bin at the LHC. Therefore, regeneration will 
dominate the observed J/ip. As a consequence, the J/ip 
Raa, V2 and (pf) at the LHC will follow the regenera- 
tion calculations in Figs. [TJ [5] and[3J 
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FIG. 3: The J/ip elliptic flow in ^sWn = 200 GeV Au + Au 
collisions at impact parameter 6 = 7.8 fm. The calculations 
are the same as in Fig. [TJ 
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FIG. 4: The ratio of regenerated to initially produced J/ip 
yields at RHIC and LHC energies in pQCD and thermal 
scenarios. At the LHC, we choose &p P = 18.9 fib [32| . Upp = 
5740 yh [33J , to = 0.3 fm and initial temperature 680 MeV. 

In summary, by constructing and solving the trans- 
port equation for charmoniums together with hydro- 
dynamic evolution of the QGP, we presented a self- 
consistent way to calculate the J/ip transverse momen- 
tum distributions. Unlike the cases at the SPS, where 
initial production is almost the only J/ip source and 
at the LHC, where the regeneration dominates the ob- 
served J/ips, both initial production and regeneration 
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are important at RHIC. The average squared trans- 
verse momentum with thermalized charm quarks fits 
the RHIC data reasonably well. In contrast to previ- 
ous calculations [12] with only sudden regeneration at 
hadronization, the dominant initial production in semi- 
central collisions and the continuous regeneration in the 
whole QGP volume lead to a rather small J/tp elliptic 
flow at RHIC. 
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